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We report a microscopic derivation of two-component Ginzburg-Landau (GL) field theory and 
the conditions of its validity in two-band superconductors. We also investigate the conditions when 
microscopically derived or phenomenological GL models fail and one should resort to a microscopic 
description. We show that besides being directly applicable at elevated temperatures, a version of a 
minimal two-component GL theory in certain cases also gives accurate description of certain aspects 
of a two-band system even substantially far from Tc. This shows that two-component GL model can 
be used for addressing a wide range of questions in multiband systems, in particular vortex physics 
and magnetic response. We also argue that single Ginzburg-Landau parameter cannot in general 
characterize magnetic response of multiband systems. 
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I. INTRODUCTION 

Ginzburg-Landau theory of single-component super- 
conductors has historically proven its strong predicting 
power and its extraordinary value as a phenomenolog- 
ical tool. This is despite the fact that formally it can 
be justified only in some cases in a very narrovir band of 
temperatures. The temperature on the one hand should 
be high enough to permit an expansion in a small order 
parameter. On the other hand the temperature should 
not be too high because the mean-field theory becomes 
invalid near Tc due to critical fluctuations. Nonetheless 
the great success of the GL theory is due to the fact that 
it yields a qualitatively correct picture in extremely wide 
range of temperatures even when its application cannot 
be justified on formal grounds. 

Shortly after the theoretical proposal of two-band su- 
perconductivity and more recently, two-component GL 
(TCGL) expansions were done in application to two- 
band systems, see e.gii. However in contrast to single- 
component GL theory, the conditions under which TCGL 
model is valid is still widely believed to be an open ques- 
tion. In this work, we resolve this question. By taking 
advantage of the recently calculated normal modes and 
length scales in two-band Eilenberger model^ we present 
the first self-consistent microscopic analysis of the ap- 
plicability of TCGL theories to describe both linear and 
nonlinear responses of two-component superconductors. 
The results validate applicability to TCGL for studying 
of a wide spectrum of physical questions, including as- 
pects of physics in low-temperatures regimes. 

The problem: the temperature range of validity of 
TCGL model is bounded from below by a requirement 
that the field amplitudes should be small. More impor- 
tant bound follows from the observed, under certain con- 
ditions, disappearance of one of two fundamental length 
scales governing the asymptotical behavior of the super- 
fluid density in microscopic theories^. This implies that 



classical two-component field theory obtained in power- 
law expansion does fail at low temperatures or at sub- 
stantially strong interband couplings^. Also, like for its 
single-component counterpart, the region of validity of 
TCGL expansion is bounded from above by a fluctua- 
tion region. 

The key difference between TCGL and single- 
component GL theory is the fact that the former has sev- 
eral coherence lengths. The existence of multiple length 
scales which arise from hybridized normal modes of the 
linearized TCGL theory, can dramatically affect the mag- 
netic response of the system^. Under certain conditions it 
results in situations where the London penetration length 
falls between two coherence lengths^, which was recently 
termed "type-1.5" regime^ (for a recent brief review 
see^). However because the condensates in the bands are 
not independently conserved, in the limit T — ^ Tc there 
should be indeed only one divergent length scale associ- 
ated with density variations. This in turn implies that 
in certain cases there also could be a temperature range 
close to Tc where long- wavelength physics is well approx- 
imated by a single-component GL theory (although this 
regime is not generic since its width is controlled by dif- 
ferent parameters than the fluctuation region and thus 
it can be non-existent because of critical fluctuations). 
Therefore the conditions of the applicability of TCGL 
are quite different from that of single-component GL the- 
ory and warrant a careful investigation, which we present 
below. 



II. MODEL AND BASIC EQUATIONS 



Expansion in powers of gradients and gap func- 
tions of microscopic equations yields the two-component 
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Ginzburg-Landau (TCGL) free energy density: 
F = (a,|A,p + ||A,|^ + X,|r»A,p) 

3 = 1-2 ^ / 

- 7(AiA; + A2AI) + |^}- (1) 

where £> = V + iA, A and B are the vector poten- 
tial and magnetic field and Ai 2 are the gap functions 
in two different bands. Despite the fields A1.2 are often 
called "two order parameters" in the literature, below 
we avoid this terminology since it is not quite accurate. 
First there is only U{1) local symmetry in this model in- 
spite of the presence of two components, since the other 
global U{1) symmetry is explicitly broken by the terms 
7 (A1A2 + A2AI). Second, and more important circum- 
stance is that the applicability of the Ginzburg-Landau 
or Gross-Pitaevskii classical field theory does not in gen- 
eral require any broken symmetries. The simplest exam- 
ple being two-dimensional supcrfluids at finite tempera- 
ture: they can be indeed be described by Gross-Pitacvskii 
classical complex field, yet they do not possess sponta- 
neously broken symmetry. Likewise in superfluid turbu- 
lence there is even no algebraic long-range order yet the 
system can be described by a classical complex field. In- 
deed in some cases, like e.g. in U{1) x U{1) superconduc- 
tors or supcrfluids one can write down two-component 
classical field theory on symmetry grounds. A C/(I) sys- 
tem such as two-band superconductors can also under 
certain conditions be described by two-component clas- 
sical field theory, although it does not automatically fol- 
lows from its symmetry. The main aim of this paper is to 
analyze under which conditions two-band superconduc- 
tors are described by TCGL theory. 

To verify applicability of TCGL theory we present 
a comparative study of linear response and non-linear 
regime in TCGL and exact microscopic theories. We con- 
sider the microscopic model of clean superconductor with 
two overlapping bands at the Fermi level^. Within quasi- 
classical approximation the band parameters character- 
izing the two different cylindrical sheets of the Fermi sur- 
face are the Fermi velocities Vpj and the partial densities 
of states (DOS) lyj, labelled by the band index j = 1, 2. 

It is convenient to normalize the energies to the critical 
temperature Tc and length to ro = hVpi/Tc- The vector 
potential is normalized by (/)o/(27rro), the current density 
normalized by c4>q/ (Sir'^r^) and therefore the magnetic 
field is measured in units (j)^) j {2i:rf^) where (/>o = nhc/e 
is the magnetic flux quantum. In these units the Eilen- 
bcrger equations for quasiclassical propagators take the 
form 

VFjnpDfj + 2Lo„.fj - 2A,5, = 0, (2) 
VFjnpD*f+ - 2c^„/+ + 2A*.9j = 0. 

Here vpj — Vfj/Vfi, w„ = {2n + 1)ttT arc Matsubara 
frequencies, the vector np = (cosflp, sinflp) parameter- 
izes the position on 2D cylindrical Fermi surfaces. The 



quasiclassical Green's functions in each band obey nor- 
malization condition + fj = 1. 

The self-consistency equation for the gaps is 

„27r 

A,; = rV / X,jfjd9p. (3) 

n=0 

The coupling matrix Xij satisfies the symmetry relations 
n-iXi2 = ?^2'^2l where n,; are the partial densities of states 
normalized so that ni + n2 = 1. The vector potential 
satisfies the Maxwell equation V x V x A = j where the 
current is 

j = -T ^ ^ /m / npgjdOp. (4) 

j=l,2 Ti=0 

The parameters aj are given by aj = A'KpnjVpj and 

p={2e/cf{roVFifvQ. 

The derivation of the TCGL functional ([IJ from the 
microscopic equationsi formally follows the standard 
scheme (we present it in the Appendix First we 
find the solutions of Eqs.(I2|) in the form of the expan- 
sion by powers of the gap functions amplitudes Ai^2 and 
their gradients. Then these solutions are substituted to 
the self-consistency equation for the gap functions which 
yields the TCGL theory. The derivation of coefficients in 
the expansion ([T|) is presented in the Appendix 1X1 Here 
we denote the values of coefficients obtained from micro- 
scopic theory as Oj^ , 6^ , A'^ and 7 which arc given by the 
expressions 

ai= pn^{\u+\T^T - Gc) (5) 
7 = pnin2Aj/DetA 
hi = pniX/T"^ 
Kt = pvF,h/4: 

where Aj = A2i/ni = Ai2/n2. Here X = 7C(3)/87r^, 

A,j- = and Gc = [TrA - VTrA^ - 4DetA]/(2DetA). 

Note that in general the derivation of the TCGL model 
is not implemented as an expansion in powers of a single 
small parameter r = (1 — T/Tc) but the outlined above 
procedure is based on the assumption of smallness of sev- 
eral parameters (gap functions and their gradients). In- 
deed the formal justification of these assumptions is not 
straightforward and to the present moment it has been 
abseniji. In the present work we show under what condi- 
tions these assumptions are rigorously justified. 

III. ASYMPTOTIC BEHAVIOUR OF THE 
FIELDS AND COHERENCE LENGTHS. 

First we investigate the asymptotical behaviour of the 
superconducting gaps formulated in terms of the linear 
modes of the density fields both in GlJ- and microscopic^ 
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theories. To find the linear modes we rewrite the equa- 
tions in terms of the deviations of the gap fields from 
their ground state values: = A^o + A^ where i = 1,2. 
To illuminate the qualitatively important physics we con- 
sider a one-dimensional case in the absence of magnetic 
field. Let us rewrite the TCGL equations by keeping on 
the left hand side the terms linear in deviations A^ while 
collecting the higher-order nonlinear terms in the r.h.s.: 

[Kid^dx^ ~ai- 36iA?o] Ai 7A2 = iVi (6) 
[K2S/dx' - 02 - 352A2o] A2 + 7A1 = 7V2. 

The r.h.s. gives nonlinear source terms Ni = 6i(3AioAf 4- 
A?). The solution of Eq.® can be found in Fourier 
representation to have the form 

K,{k) = R-'N,{k) (7) 

where 

^12 = -R21 = 7 
R^, = - [K,e + a, + 3fa,Afo] . 

In this case the response function has two poles in 
the upper complex half-plane k = i/i/f and k = i/i^ which 
determine the two inverse length scales or, cquivalently, 
the two masses of composite gap functions fields^, which 
we denote as "heavy" and "light" (i.e. hh > Ml)- 
Let us set Ki = K2 which can be accomplished by 
rescaling the fields Ai 2. Then the matrix R~^{k) can be 
diagonailized with the fc-independent rotation introduc- 
ing the new linear modes of the fields $/3 = UpiAi and 
the sources = UpiNi where /3 ~ L,H and i = 1,2. 
The rotation matrix U is characterized by the mixing 
anglei^ as follows: 

y~ sm Oh cos Oh J 

Using the diagonal form of the response function 
R^^{k) in the real-space domain we obtain 

(9) 

where = [Nfi{x) +2Np{0)]e-t'f"=dx/2nfi is cho- 
sen so that to satisfy the boundary condition $/j(0) = 
Np{0)/^Jj which corresponds to the condition Ai.2(0) = 
at a: = 0. 



A. The limit r ^ 0. 

The expression ([9]) shows that two fields ^l.h vary at 
distinct coherence lengths: S^h = l/^J■H and = I/a^l- 
They constitute fundamental length scales of the TCGL 
theory ([T]). They characterize the asymptotical relax- 
ation of the linear combinations of the fields Ai 2, the lin- 
ear combinations are represented by the composite fields 



^L,H- Our calculation shows that these length scales 
behave qualitatively different in the limit r — >■ 0. In- 
finitesimally close to Tc the largest length diverges as 
S^L ^ T~^/^ while the smaller ^h remains finite. Sim- 
ilar behavior also follows from full microscopic calcula- 
tion shown on Fig.([T])b,c,d where the temperature de- 
pendence of masses fiL.H is plotted. The presence of the 
non-diverging length scale makes the qualitative dif- 
ference with the single-band GL theory but indeed does 
not contradict the standard textbook picture that in the 
limit r — > the mean field theory of a J7(l) system should 
be well approximated by single-component GL model. 
As we show below the amplitude of the "heavy" mode 
vanishes in the t ^ Q limit faster than that of "light" 
mode. Neglecting the "heavy" mode contribution one 
indeed obtains a single-component GL theory. 

We can use the Eq.® to evaluate the asymptotical 
amplitudes of ^h,l(j') in terms of the powers of the ex- 
pansion parameter r, in the limit r — >■ 0. The goal is 
to evaluate how the contributions from different length 
scales affect overall profile of the fields as they recover 
their ground state value away from x = 0. First we note 
that the source terms Nl^h{x) are confined at the region 
determined by the coherence length x < ^l- Inside this 
region the amplitude of the deviations of the gaps from 
the ground state values are large so that 

A,(x)^Ao.^t1/2. 

Thus the amplitude of sources is of the order Nl^h ^ 
T^l"^ . Let us consider the first term in the expression ([9]) 
for /3 = L at the asymptotical region x > In this case 
the integration is confined within xi < and yields the 
following estimate $l ~ A^e"'^'^^ where 

A c ^3/2 ^1/2 

The second term in the Eq.([5]) gives the contribution of 
the same order to the amplitude of the "light" mode. 

The amplitude of the "heavy" mode is determined en- 
tirely by the first term in Eq.® for /3 = H. We consider 
the asymptotical region a; therefore in this estimate 

we put Nh{xi) ~ Nh{0) so that 

<fH{x) = NH{0) {2e-^"-^l) /njj. 

The function ^h (x) has a characteristic scale £,h ~ 1/ M// 
and its overall amplitude is determined by the factor 

Ah - Nh{Q)/A - r3/2. 

Thus in the limit r — > the "heavy" mode drops out be- 
cause of the vanishing amplitude Ah ~ t^^"^ as compared 
to the "light" mode Al ^ t^I"^ . Note that it is thus prin- 
cipally incorrect to attribute different exponents directly 
to the functions A^ and to assume that they become equal 
in the limit t — ^ as claimed in^ and followed in some 
other literature^. 

On the qualitative level we give a less technical but on 
the other hand more intuitively transparent description 
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of the limiting behavior of the fields near T^.. We consider 
the situation where the coefficients of quadratic terms in 
Eq. (IlJ can be written in the form Oj (T) = aj [T—Tj ) with 
aj > 0. Thus in the small t limit, first the weakest super- 
conducting component becomes passive: it has nonzero 
superfluid density only because of the bilinear Joseph- 
son coupling 7(AiA2 +c.c.). To elucidate what happens 
in the r — > limit one can redefine fields using the fol- 
lowing transformation A = XD4'_D(r) +Xs5's(r) where 
Xd = (q!2(Td — T2), 7) , Td = Tc is the critical temper- 
ature, Xs = (7, Qfi {Ts - Ti)f and Ts = Ti+T2-To < 
To- This transformation, mixes the gaps and produces 
a representation where the bilinear Joscphson coupling 
between the new fields is eliminated at the cost of intro- 
ducing mixed gradient and fourth-order couplings: 

T.,,,Ki=D,s ft,fe/i?e(*.*,vi/*vi/;). (10) 

In the limit r — > the dominant component ^E'o be- 
comes single-component GL order parameter while the 
component ^'s is passive. However instead of being in- 
duced by the bilinear Josephson coupling, s is induced 
by the terms like *s*dI*oP- Thus for r -)■ one has 
"Hd ^ T^/^ From the fact that /3s (T - Ts) and Psddd 
arc finite at r ~ it follows that in the same limit we 
have D ^ ^ s ^ t^^"^ . It means that in the limit r — )• 
retaining only the terms containing ^'^i is justified and it 
allows to approximate the model by a conventional single 
order parameter theory. 

For infinitesimally small r there remains only single 
GL equation for the order parameter 

(11) 

where Pdddd = {alo^i + ain^2)/ai and 

13d = > (12) 

_ aiDK2 + a2DKi 
aia2{iD - J-2) 

where ai,2D = ai^2{TD)- 

Thus it is the disappearance of the amplitude of the 
subdominant mode which allows one to take a single- 
component GL limit in this mean-field theory. The 
"heavy" mode with finite mass even infinitesimally close 
to Tc is generated here by the presence of the term 
l-D^'sP ~ Cff^l^sP where S.H is a finite length not di- 
verging at Tc- Therefore near Tc one has jD^I'sp ^ 
^H^l^-sp ^ ''"'^ which is of the same order as the other 
terms in the free energy functional. Note that the fields 
^'_D,s introduced here are not directly related to the nor- 
mal modes of the system since the mixed gradient and 
quartic terms will lead to mode mixing^, at finite r. Also 
note that despite there is a growing disparity of the co- 
herence lengths at small but finite r when one approaches 
critical temperature, it does not imply that one necessary 



falls into type- 1.5 regime because < A < ^2 is only nec- 
essary but not sufficient condition for the appearance of 
this regime. That is, a system in some of these cases is 
type-I despite having ^1 < A < ^2- 



B. GL theory at finite r. 

Unfortunately the limiting r — > analysis docs not 
have much physical significance in a generic two-band 
system. First, the mean-field theory becomes invalid 
in the same limit t — > so the regime where the sys- 
tem is well described by single-component GL theory 
can be cutoff by critical fiuctuations. More importantly, 
as we argue below, this analysis is in general inapplica- 
ble for an assessment of, e.g. magnetic response of the 
system. The magnetic response is a finite-length scale 
property and requires finite-r theory. Finally, as shown 
in microscopic calculations the masses of the fields in 
two-band models in certain cases change rapidly and in 
a non-trivial way with decreasing temperature. Thus 
a limiting t — ^ analysis in general cannot give even 
an approximate physical picture even at very small r. 
In particular that implies that in a two-band system, a 
Ginzburg-Landau parameter (which one, may in princi- 
ple construct in the r — > limit at a mean-field level) 
is not a useful characteristic. Rather it is required to 
make an accurate quantitative study of two-band the- 
ory at finite r to determine the conditions under which 
the model can be described by singe- or two- compo- 
nent GL theory or does not allow a description by any 
such GL functionals at all. In order to do it we uti- 
lize the exact form of the response function (i.e. 
valid at any T) found from the linearized microscopic 
theory according to the procedure developed in^. In con- 
trast to the GL theory, the microscopic response func- 
tion R^^(k) has branch cuts along the imaginary axis 

starting at point zkikt,c where k^c = 2 yA^j+T^'^^ ("wc 
assume that A02 < Aqi). Inside the circle |fc| < \k},c\ 
shown by the white area in Fig.([T])a the response func- 
tion is meromorphic, i.e. its singular points are only poles 
shown by the red crosses. The non-meromorphic region 
is marked by the yellow shade in all panels of Fig.([T]). In 
general inside the meromorphic circle there can be two 
poles of R~^{k) at Im{k) > shown by red crosses in 
Fig.(IT|)a. Analogously to TCGL these poles determine 
the masses of the "heavy" and "light" modes, and 

thus the corresponding coherence lengths. The contribu- 
tion of the branch cut contains the continuous spectrum 
of length scales shorter than 1 / k^c which can not be de- 
scribed within GL theory. Moreover for some parame- 
ters (e.g. at strong Josephson coupling) one of the poles 
which corresponds to the "heavy" mode can disappear 
by merging with the branch cut. In this case there is 
only one fundamental length scale left since the contri- 
bution of the "heavy" mode can not be separated from 
the branch cut. 

The microscopically calculated temperature dependen- 
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cies of masses of the modes in a superconductor with 
weak interband couphng are shown in Fig.(IT])b by red 
sohd hues. For a reference we also plot masses in 
C/(l) X C/(l) theory which has two independently diverg- 
ing coherence lengths at T = Td and T = Tc2 (chosen 
to be Tc2 = O.STci). For coupled bands the hybridiza- 
tion of modes removes the divergence at T = Tc2 and 
introduces the avoided crossing point of the "heavy" and 
"light" modes. 

Let us now assess the applicability of minimal TCGL 
model Eq.(IT|) without using expansion in powers of t. 
Compared to the previous worksi, we use more com- 
plicated temperature dependence of the coefficients de- 
rived in the Appendix [X] Let us compare the behavior 
of the masses of the modes in the microscopically de- 
rived TCGL and a full microscopic theory. It is shown 
for the cases of weak and strong interband coupling in 
Fig.dDlCjd. We have found that TCGL theory describes 
the lowest characteristic mass /xl(T) with a very good 
accuracy near T,. [compare the blue and red curves in 
Fig. ([T|)c,d]. Remarkably, when interband coupling is 
relatively weak [Fig.([T])c] the "light" mode is quite well 
described by TCGL also at low temperatures down to 
T = 0.5Tc around which the weak band crosses over from 
active to passive (proximity-induced) superconductivity. 
Indeed the r parameter is large in that case and cannot 
be used at all to justify a GL expansion. Nonetheless 
if the interband coupling is small one does have a small 
parameter to implement a GL expansion for one of the 
components. Namely one can still expand, e.g. in the 
powers of the weak gap | A2 I/ttT ^ 1. On the other hand 
for the "heavy" mode we obtain some discrepancies even 
relatively close to Tc, although TCGL theory gives quali- 
tatively correct picture for this mode when the interband 
coupling is not too strong. More substantial discrepan- 
cies between TCGL and microscopic theories appear only 
at lower temperatures or at stronger interband coupling 
[Fig.([T])d] where the microscopic response function has 
only one pole, while TCGL theory generically has two 
poles. 

Comparison of the masses of normal modes of a ?7(1) x 
C/(l) (dotted blue lines) and a weakly coupled two-band 
U{1) model (red lines) shown in Fig.(IT])d demonstrates 
that adding a Josephson coupling removes divergence of 
coherence length of the weak band. This is because the 
Josephson coupling represents explicit symmetry break- 
down from U (1) X U{1) to U{1) and thus eliminates one of 
the superconducting phase transitions at lower Tc- How- 
ever when this coupling is weak one of the coherence 
lengths has a substantial peak around that temperature. 
The peaked behaviour of coherence length near the criti- 
cal temperature of the weak superconducting band is has 
clear physical manifestation in the temperature depen- 
dence of the vortex core size. Let us note that to assess 
the overall size of core requires analysis of full nonlinear 
theory. In Fig.© we plot the sizes of the vortex cores 
in weak and strong bands calculated in the full nonlinear 
model according to the two alternative definitions. The 



first one is the slope of the gap function distribution at 
r = which characterizes the width of the vortex core 
near the center R^j = (din Aj/dr)"i (r = 0) [Fig. ©a]. 
The second one is the healing length Lhj defined as 
Aj{Lhj) = 0.95Aoj [Fig.Qb] (i.e. this length is not di- 
rectly related to exponents but quantifies at what length 
scales the gap functions almost recover their ground state 
values). Both definitions demonstrate the stretching of 
the vortex core in the weak component related to the 
peak of the coherence length shown in the Fig.([T])d. Note 
that the weak band heahng length Lh2{T) in Fig.Qb 
has maximum at the temperature slightly larger than 
Tc2 which is consistent with the fact that the maximum 
of coherence length (equivalently the minimum of the 
field mass /i^) in Fig.(IT|)d is shifted to the temperature 
above Tc2 (Tc2 is defined as the lower critical temperature 
in the limit of no Josephson coupling). 



C. Effects of higher order gradient terms. 

The origin of the small disagreement between the 
TCGL and microscopic masses of the "heavy" mode is 
the absence of higher-order gradient terms in the expan- 
sion IH) . The inclusion of higher order gradients means 
adding more terms to the Taylor expansion of the func- 
tion R{k) which is known to converge inside the circle 
|fc| < |fcf,c|. Using this procedure one can get a bet- 
ter agreement with microscopic theory (compare green 
dashed and red sohd lines in Fig.(IT|)c,d). However, such 
an extension is hardly useful because as a byproduct 
it generates unphysical artifacts such as many parasitic 
poles of the response function R~^{k) outside the circle 
\k\ < \khc\- The position of parasitic poles appearing in 
the sixth-order gradient expansion is shown by green cir- 
cles in the Fig.([T^). These parasitic poles lie outside the 
imaginary axis, thus yielding unphysical oscillating con- 
tributions to asymptotical behavior of the corresponding 
linear modes. 



D. Characteristic length scale of the phase 
difference variations. 

In the U{1) X U{1) system one has two massive modes 
associated with the modules of the complex fields and 
a Goldstone boson associated with the phase difference. 
If one adds a Josephson coupling there appears also the 
third mass parameter. When Josephson term is present 
the phase difference acquires a preferred value. Its de- 
viations from the preferred value are characterized by a 
mass parameter. In the constant density approximation 
the terms of the TCGL functional which describe the 
phase difference mode are: 

I KA^'t'fh (^(^^ - S2)r 7A10A2O COS(01 - 62) 

(13) 
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Expanding the second term gives the mass parameter for 
the phase difference mode 



^{01-02 
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(14) 



It is useful to consider the behavior of this mass in the 
hmit T — )- Tc- In that case we have Aio,20 oc i/l 
Thus 



hm (m(ej_e2)) — const 



(15) 



To summarize this part we have shown that the TCGL 
model of the form given by Eq. ([T|) with microscopically 
derived temperature dependencies of coefficients is over- 
all highly accurate at elevated temperatures. The small 
discrepancies with microscopic theory affect only short- 
length scale physics which implies that TCGL model 
gives the precise answer for long-range intervortex forces. 
Also we find that in some cases the TCGL model provides 
an accurate description of the large length scales physics 
at temperatures much lower than Tc. In Appendix lA II 
we discuss the origin of the disagreements between these 
results and some of the recent literature^. 



IV. VORTEX STRUCTURE: TCGL VS 
MICROSCOPIC THEORY. 

For inhomogcneous situations, such as vortex solu- 
tions, the overall profiles of the fields is affected not only 
by fundamental length scales (i.e. coherence lengths) but 
also by nonlinear effects. 

Let us now study non-linear effects case of vortex so- 
lutions. Obviously, because of the growing importance of 
nonlinear effects at lower temperatures the Eq.([T]) cannot 
describe quantitatively well the total structure of vortices 
when T <C Tc- In Fig.(l3|) we compare the vortex solu- 
tions in the self-consistent microscopic theory (red dot- 
ted curves) and in the corresponding TCGL theory with 
coefficients obtained by expansion (blue dashed-dotted 
curves). One can sec that at elevated temperatures the 
agreement is very good but for lower temperatures there 
is a growing discrepancy. One of the reasons behind the 
discrepancy is the trivial shift of the ground state values 
of the fields by nonlinearities. Note that at the level of 
GL theory the inclusion of more nonlinear terms merely 
renormalizes masses and length scales but does not alter 
the form of linear theoryi. Thus in the current exam- 
ple of the full nonlinear model it is also reasonable to 
check if one could get a better agreement with micro- 
scopic theory by treating the coefficients in the minimal 
TCGL model Eq.([T]) phenomenologically. For all practi- 
cal purposes this provides alternative route to the more 
restrictive approach of finding a refined microscopic ex- 
pansion. A good agreement with the microscopic theory 
in this procedure will imply that the system does posses 
a description in terms of a classical two-component field 
theory. 



T 


0.98 


0.8 


0.7 


0.5 


0.4 


0.2 


bi 


0.95bi 


0.85fei 


0.765i 


0.545i 


0.446i 


0.18bi 




Ki 


0.8i?i 


0.657?i 


0.5A'i 


0.4Ai 


O.lSA'i 


b2 


b2 


b2 


b2 


0.7662 


0.6462 


0.2962 


K2 


K2 


0.55i?2 


0.357?2 


0.3A'2 


0.15A2 


O.O8A2 



TABLE I: Fitting of TCGL coefficients to match the solutions 
of exact microscopic equations. Here we denote the values of 
coefficients obtained from microscopic theory as a^, bv, 
and 7. 



We compared the vortex solutions in the TCGL theory 
with fitted coefficients and the exact microscopic model 
for the particular example of the system with coupling 
constants An = 0.5, A22 = 0.46, A12 = A21 = 0.005 and 
VF2/VF1 = 5. The values of the coefficients which provide 
the best fit are listed in the table ID 

By the green dashed lines we show the fits obtained by 
setting the values of the TCGL coefficients as listed in the 
table m By doing it we find that the solutions of micro- 
scopic equations in a large region of parameters can be fit- 
ted with excellent accuracy by the effective TCGL theory, 
even at quite low temperatures. The main discrepancies 
at very low temperatures arise due to non-local effects 
which lead to the disappearance of the "heavy" asymp- 
totic mode as well as due to Kramer-Pesch-like vortex 
core shrinking^iiS in both components which cannot be 
captured in the TCGL field theory. Note however that 
even in the case where TCGL description starts break- 
ing down, the discrepancy is mostly pronounced near the 
origin of the core, while the soft modes and long-range 
intervortex interaction can be well described by a phe- 
nomcnological TCGL theory. 



V. CONCLUSIONS 

TCGL model is widely used for describing various 
aspects of multiband superconductivity. However the 
TCGL expansion has never been rigorously justified for 
two-band systems, and the current literature contains di- 
ametrically opposite claims regarding the validity of the 
expansion or such basic aspects as the form of TCGL 
functional and behavior of the coherence lengths near 
Tjii^. We investigated under which conditions a two- 
band system can be described by a TCGL theory. First 
we obtained a TCGL model with a microscopically de- 
rived temperature dependence of coefficients (more gen- 
eral than what could be obtained in a straightforward 
T expansion) and demonstrated that it gives an accu- 
rate description of length scales and vortex solutions at 
elevated temperatures by a comparison with an exact mi- 
croscopic theory. Second we have shown that, in a much 
wider range of temperatures, the minimal TCGL model 
with phenomenologically adjusted coefficients gives an 
accurate description of linear and nonlinear physics such 
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as vortex excitations and thus the magnetic response of 
the system. 

The existence of two coherence lengths , £,h along 
with the magnetic field penetration lengths A in the 
TCGL model makes it impossible in general to define the 
Ginzburg-Landau parameter in two-band systems, unless 
one takes the limit r — >■ 0. In contrast to U{1) x U{1) 
superconductors, the two-band systems have only U{1) 
symmetry, and as we discussed above it guarantees that 
one of the modes drops out of the mean field theory 
the limit r ^ allowing one to define in that limit 
KGL = idT Tc)/\{T -> Tc). However, even slightly 
away from the limit r — ^ 0, when interband coupling is 
weak the ratio / A has a very strong temperature depen- 
dence and the second mode develops with the coherence 
(^H- Thus in general kql cannot be used as universal 
characteristic of magnetic response of two-band systems. 
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Sciences, Swedish Research Council, US NSF CAREER 
Award No. DMR-0955902, "Dynasty" foundation. Pres- 
idential RSS Council (Grant No. MK-4211. 2011.2) and 
Russian Foundation for Basic Research. 



Appendix A: Microscopic model and derivation of 
TCGL 



(A5) 



Expressing GLi from the equations above we obtain 

n,GL, =n,{^-G]^,- ^^A, (A6) 



VDetA 



DetA 



n2GL2 = 112 ( - G ) A2 - ^^^Ai (A7) 

VDetA / DetA 

Comparing these Eqs with Eqs ([5]) we obtain the expres- 
sion for the coefficients 

a,= pn,{\,+\nT -Gc) (A8) 
7 = (071171,2 A j/Det A 

where Aj = A2i/7t,i = Xi2/n2- The temperature is nor- 

Kj 

TrA - VTrA2 - 4DetA 



mahzed to the T^. Here X = 7C(3)/87r2, A^ A,^-^ and 



Gc = G{T, 



2DetA 



We have used the expression G(T) = G{Tc) - InT. Near 
the critical temperature InT « — r and we obtain = 
?i,Aj(T - Ti) where T, = (1 + G^ - A,,)- 



1. Ginzburg-Landau expansion. 

To derive differential GL equations ^ from the mi- 
croscopic theory, first we find the solutions of Eilenberger 
Eqs.® in the form of the expansion by the gap functions 
amplitudes |Ai_2| and their gradients |(iI>np)Ai_2|. Then 
these solutions are substituted to the self-consistency 
Eq.dSj. Using this procedure we find the solutions of 
Eqs. ([2]) in the form: 



/ 



|ApA 
2ul 



^(r>np)A + ^(Z)np)(£)np)A. 



2uj\ 



(Al) 

and /+(np) = /*(— Hp). Note that a GL expansion is 
based on neglecting the higher-order terms in powers of 
|A| and |(£>np)A|. Indeed this approximation naturally 
fails in in a number of cases. In this work we determine 
the regimes when it can be justified, in particular by di- 
rect comparison with exact microscopic model. Let us 
determine microscopic coefficients in the GL expansion. 
Substituting to the self-consistency Eqs.Q and integrat- 
ing by Op we obtain 

Al = (AiiAi + Ai2A2)G + (AiiGLi + Ai2Gi2) (A2) 

A2 = (A21A1 + A22A2)G + (A21GL1 + A22Gi2) (A3) 



where 



G 



(A4) 



2. Remarlt on r-expansion 

Here we comment on the origin of the qualitative dis- 
agreement of our result compared to the worki which 
aims at calculating higher order corrections in r = (1 — 
T/Tc). Fist let us make a few general remarks: Note that 
in the derivation of TCGL theory we do not implement 
an expansion in powers of r = (1 — T/Tc). Instead we 
retain more complicated temperature dependence of the 
coefficients. Also we stress that any approach to GL ex- 
pansion depends on what parameters are assumed to be 
small, the question is always how and for what param- 
eters such an assumption is justified. The origin of the 
principal difference in the behavior of the length scales 
(Refi^ asserts that there are two divergent length scales 
when r — > 0). It originates in the adoption in^ of a 
C/(l) X U{1) theory as the leading order in expansion fol- 
lowing the erroneous derivation in^ (see the discussion 
of the errors in that derivation in^). Another problem 
with a straightforward implementation of the expansion 
by r is that in general it is incontroUable in the next to 
leading order in two-band theories, if one explicitly re- 
tains two gap fields. This is because in contrast to the 
single-component GL theory, in general it is not possible 
to classify different terms by powers of the parameter r. 
As shown in the main body of the paper, system contains 
a mode with non-diverging coherence length so that the 
spatial derivatives in general do not necessary add the 
power of T. Also since the worki uses as a leading order 
the incorrect derivation from^ it requires adjustments. 
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Appendix B: Asymptotic in two-dimensional vortex 
problem 



The consideration of asymptotic modes in Sec lIIII can 
be generalized for the two-dimensional axially symmet- 
ric problem, which allows to treat the asymptotical be- 
haviour of the gap functions far from the vortex core. 
First of all in this case one should substitute the <P /dx^ 



by = d^/dr^ 



^d/dr in Eq. ©. 



Choosing the proper value of the mixing angle the l.h.s. 
of Eq. ([5]) can be diagonalized and the system acquires 
the form 



(Bl) 



where the nonlinear Nh.l are obtained according to the 
rule dl]). 

Our interest is the asymptotical behaviour of the fields 
^L,H determined by the equation above. The solution of 
Ea. (|Bip can be found in Fourier representation A(fc) = 
K{x)e'^^^ dx to have the form ([7]). In this particular 
case the response function is a diagonal matrix: 

In the real-space domain the field components can be 



expressed with the help of Fouricr-Besscl transform 



= / Mkri)jQikr)R^^Nj{ri)kdkridri. 
Jo 

The integration by k in this expression can be performed 
by transforming the contour in the complex plane. Using 
the exact form of the response function the fields asymp- 
totic is found to be given by the following expression 



■^^{r)^^TKo{^l^r) ridriIo{ii,ri)N,{ri) + (B2) 



7r/o(^ir) / ridriKQ{^iri)Ni{ri) 



where -R'o and /q are modified Besscl functions having 
the following asymptotics A'o,/o(2;) ~ e^^ j^/x. 

The expression (jB2|) yields a number of length scales 
characterizing the asymptotical relaxation of the gap 
fields. The largest length is the mean field coherence 
length = 1//^L ~ However the presence of the 

another linear mode in the theory sets the scale which 
is proportional to S^n = l/fiH- This scale remains finite 
even at T = Tc but its amplitude vanishes. 
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FIG. 1: (a) Comparison of the response function singulari- 
ties in the complex k plane given by the exact microscopic 
and microscopically derived TCGL theories. Red crosses are 
the physical poles of microscopic theory. Blue squares corre- 
spond to the conventional TCGL theory while green circles 
show the parasitic poles appearing in TCGL expansion up to 
sixth order in gradients. The white circle is the area where 
R{k) is analytical and R~^{k) is meromorphic. In all plots 
(a-d) the yellow shade indicates the area where the response 
function is not meromorphic. (b) Comparison of the masses 
of normal modes of a [/(I) x (7(1) (dotted blue lines) and 
a weakly coupled two-band U{1) model (red lines). For the 
weak component in the (7(1) x (7(1) we also plot the corre- 
lation length for superconducting fluctuations in the normal 
state for T > Tc2- As clearly seen in that plot, adding a 
Josephson coupling removes divergence of coherence length of 
the weak band. This is because the Josephson coupling repre- 
sents explicit symmetry breakdown from (7(1) x (7(1) to (7(1) 
and thus eliminates the phase transition at lower Tc. However 
when this coupling is weak one of the coherence lengths has 
a peak around that temperature, (c) and (d) Comparison of 
field masses given by microscopic (solid lines), TCGL (dotted) 
and TCGL with sixth-order gradients (dashed lines) theories. 
The microscopic parameters are An — 0.5, A22 = 0.426 and 
A12 = A21 — 0.005; 0.01; 0.1 for (b,c,d) correspondingly, (e-f) 
Comparison of the mixing angle behaviour given by the ex- 
act microscopic (red lines) and microscopically derived TCGL 
theories (blue line). Parameters are the same as on the pan- 
els (c-d) correspondingly. Panels (d) and (f) show a pattern 
how the TCGL theory starts to deviate from the microscopic 
theory at lower temperature when interband coupling is in- 
creased. 




FIG. 2: (a) Sizes of the vortex cores i?ci,2 and (b) healing 
lengths Lhi,2 in weak (blue curve, open circles) and strong 
bands (red curve, crosses) as functions of temperature. The 
parameters are An — 0.5, A22 ~ 0.426, A12 — X21 = 0.0025 
and VF2/VF1 ~ 1. In the low temperature domain, the vor- 
tex core size in the weak component grows and reaches a 
local maximum near the temperature Tc2 (the temperature 
near which the weaker band crosses over from being active to 
having superconductivity induced by an interband proximity 
effect)^. In the absence of interband coupling there is a gen- 
uine second superconducting phase transition at Tc2 = 0.5rci 
where the size of the second core diverges. When interband 
coupling is present it gives an upper bound to the core size in 
this temperature domain, nonetheless this regime is especially 
favorable for appearance of type-1.5 superconductivity^. 
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FIG. 3: Behavior of the gap functions Ai,2('') in a vortex 
solution. Comparison of the results of exact microscopic cal- 
culation (red dotted lines), TCGL with microscopically cal- 
culated coefficients (blue dash-dotted lines) and TCGL with 
phenomenologically fitted coefficients (green dashed lines) at 
(a) T = 0.98, (b) T = 0.8, (c) T = 0.7, (d) T = 0.5, (e) 
T = 0.4, (f) T = 0.2. Couphng constants are An = 0.5, 
A22 = 0.46, A12 = A21 = 0.005 and vf2/vfi = 5. 



